Mode Entanglement and Entangling power in Bosonic Graphs 
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We analyze the quantum entanglement properties of bosonic particles hopping over graph structures. Mode- 
entanglement of a graph vertex with respect the rest of the graph is generated, starting from a product state, by 
turning on for a finite time a tunneling along the graph edges. The maximum achieved during the dynamical 
evolution by this bi-partite entanglement characterizes the entangling power of a given hopping hamiltonian. 
We studied this entangling power as a function of the self-interaction parameters i.e., non-linearities, for all the 
graphs up to four vertices and for two different natural choices of the initial state. The role of graph topology 
and self-interaction strengths in optimizing entanglement generation is extensively studied by means of exact 
numerical simulations and by perturbative calculations 
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I. INTRODUCTION 

Generation of quantum entanglement is an essential prereq- 
uisite for the majority of quantum information and computa- 
tion tasks 1 1]. To this aim many -body interactions have to be 
cleverly engineered and carefully controlled in the variety of 
physical systems that pretend to represent viable candidates 
for quantum information processors. The resulting quantum 
dynamics can be then characterized by means of its entangling 
capabilities in different ways QIHQIsl]. These entangling 
power measures provide a quantitive to assess the dynamical 
strength |5] of a quantum evolution as physical resource for 
quantum information processing (QIP) protocols. In this re- 
spect the optimization of the entangling power in the set of 
physically available interactions clearly represents a quite rel- 
evant goal. 

In this paper we will address this problem in the context 
of bosonic systems on lattice structures e.g., ultracold bosonic 
atoms in optical lattices. This kind of systems have been al- 
ready rather extensively studied in the QIP related literature 
BlililiSlilllillilllililin view of their experimental 
feasibility and inherent great flexibility that make them almost 
ideal candidates for studing and simulating L16J a plethora 
of important quantum many-body phenomena e.g., quantum 
phase transitions iHtIi . 

In the following we will adopt the view of quantum entan- 
glement in systems of indistinguishable particles advocated 
in Refs |18, 19, 20, 21, 22, 231 in alternative to the comple- 
mentary one pursued in Refs. 124, 25, 26 1. In the former the 
subsystems that get entangled are provided by bosonic modes 
whereas in the latter are particles themselves. The mode- 
entanglement approach is an intrinsecally second-quantized 
one: the tensor product structure of the state-space necessary 
in order to define entanglement is provided the natural iden- 
tification of the bosonic (fermionic) Fock space with a set of 
linear oscillators (qubits) associated with the single-particle 
wavefunctions. 

We shall analyze systems of bosonic particles hopping 
over inhomogneous spatial structures 11271 12811 mathematically 
modelled by graphs [6]. The main objective is to determine 



the role that the graph topology plays in the entangling power 
of the quantum dynamics described by a Bose-Hubbard model 
ll29ll . For a given tunneling graph structure we will study how 
the entangling capabilities are affected by the amount of non- 
linearity i.e., bosons self-interactions present in the Hamilto- 
nian. Numerical simulations as well as perturbative analytical 
arguments will be provided for all the graphs up to four ver- 
tices. This analysis allows to optimzes the entanglement pro- 
duction, from an initial product state, obtainable by switching 
on-and off the tunneling between the various garph vertices 
i.e., one-body operator control, and by varying e.g., by Fes- 
dbach resonances, the strength of the boson-boson local cou- 
pling. We believe that the theoretical framework developed in 
this paper will have direct experimental relevance. 

The paper is organized as follows: in Sect II is briefly re- 
called the Bose-Hubbard Hamiltonian and in Sect III it will 
be applied to the case of hopping over spatially inhomogeneus 
structures (graphs). In Sect IV we will describe the notion of 
bosonic mode-entanglement for the systems of interest and in- 
troduce the functional for measuring the entangling power of 
a quantum evolution. Sect V is the major one, it contains the 
results of numerical simulations for all the graphs up to four 
vertices. Finally Sect VI provides a summary and conclusions. 



II. THE BOSE-HUBBARD MODEL 

The systems we will study in this paper are all given by L 
bosonic sites filled with N = L bosonic particles. Some of 
the sites are linked to each other in the sense that the particles 
are allowed to tunnel from one site to another The links ex- 
isting between the sites will be described later (section III) by 
means of graphs. The particles in a site interact via a two-body 
interaction process (self-interaction). 

The proper framework that can be used to describe the sys- 
tems under study is given by the Bose-Hubbard (BH) model 
ll29ll . This model has been recently use to investigate many 
interesting properties of different systems going from aiTays 
of BEC where the sites form an infinite (periodic) lattice 
10], to systems in which the number of sites is finite (dimer, 
trimer)|l7'|. Moreover this model has recently been used to 
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describe new possible schemes for quantum computing (T^. 
According to this model the Hamiltonian of the system of L 
sites and N particles can be written as: 

i i <i,j> 

where: the indexes i,j run in V — {0, 1, 2..., L — 1}, the set of 
sites; cj and Ci are bosonic creation and annihilation operators 
that create/destroy a particle in the locale mode (site) i; rii — 
cjci are the corresponding occupation number operators. 

The parameters 7^ can be used to model the offsets in 
the ground state energies of the different sites. We suppose 
7i = 7 = const \/i and from now on we neglect this term. In 
fact, since the total number of particles in our systems will be 
held fixed, the term 7 ^ • rii — -fN commutes with the Hamil- 
tonian (1) and its contribution in the dynamics of the system 
consists only by a global phase. The parameters control the 
non-linear two-body interaction between the particles in each 
site i, while the parameters Tij control the tunneling processes 
between the sites i, j. We consider only tunneling between 
next-neighbourghring sites and we suppose to be able to set 
the tunnel parameters Tij and the self-interaction parameters 
Ei to some fixed desired values. These will be held fixed dur- 
ing the evolution, i.e. we will be intersted only on the free 
evolution of our systems. The latter will be described in detail 
in the next section. 



III. GRAPHS AND HAMILTONIANS 

As described in the previous section, each of the systems 
we consider is given by a lattice of L vertices (sites or modes) 
filled with N = L bosonic particles and we suppose that 
its Hamiltonian can be described in the framework of the 
Bose Hubbard model. We consider systems in which the self- 
interaction terms are all equal {Si ~ e Vi), while for the hop- 
ping parameters we consider the case G {0, t} Vi, j. We 
focus on the free evolution of the simpliest of these systems, 
i.e. N = Le {3,4}. 

For L = = 3, 4 we consider only inequivalent systems, 
in the sense specified below, each characterized by a different 
rooted graph T = {V, E, r), where V = {0, 1,2,.., L- 1} 
is the set of vertices, E = : i,j G V, i 7^ j}is 

the set of edges and r G V is the root of the graph. The 
latter will be considered as the mode (vertex) which is rele- 
vant for the calculation of the mode entanglement. For each 
rooted graph F the pair {V, E) identifies a particular Bose 
Hubbard hamiltonian H^y^^y since the pair (V, E) identifies 
a set of sites and the way in which these sites are connected. 
With this Hamiltionian one can build the evolution operator 
'U(v,E){'t) = exp(-iii?(v,E)). 

If the pairs {V, E) and {V, E*) are such that E ^ E* they 
give rise to two sets of inequivalent rooted graphs since the 
connectivity of the two graphs they identify is different. Mo- 
rover, for a given pair {V, E) one can define a certain rooted 
graph by specifying the root site r G V, but only some of 
the possible choices can be considered inequivalent. If two 



a) 

J 


K 

J. ^ 

A 


I 


A 

\ 


L 


1^ = 1 


= z 


1 


2 


1 2 


b) 




r 











I 

1 


2 


1 


2 



FIG. 1: a) inequivalent rooted graphs for N = L = 3\ b) discon- 
nected graphs for N = L — 3. 



rooted graphs differ only for the root, they will be considered 
inequivalent if there is a difference in how the roots are di- 
rectly connected to the other vertices of the graph or if there is 
a difference in the connectivity of the subgraph F' = {V' , E') 
constituted by the set of vertices V' ~ V\r and the set of 
edges E' = & V',i^j}. 

We now give an example of the procedure used to identify 
the different inequivalent systems by starting with N = L = 
3. In this case we will consider three different inequivalent 
rooted graphs (see Fig. la). 

To build these graphs we first start by labeling the root ver- 
tex with 0, the left (L) vertex with 1 and the right (R) vertex 
with 2. The choice of the labels for the vertices of the sub- 
graph F' = {V' , E'), where V' = {L, R} is irrelevant for our 
analysis (we could for example choose L = 2 and R = 1) 
since from the point of view of the entanglement of the root 
mode they are interchangeable and the relevant feature is how 
r is directly connected to the vertices of the subgraph F' and 
the connectivity of the latter. So we can say that, for a given 
i?',the two labeling (0, 1, 2) and (0, 2, 1) are equivalent and 
consider only the case (r, L, R) = (0, 1, 2). 

The rooted graphs Fq and F2 give rise to the same "kind" 
of BH Hamiltonians, in the sense that they both have two hop- 
ping terms; the only thing that changes is the labeling of the 
modes: 

Hr„ = e{nl + nj + nl) + r(4ci + 4^2 + {h.c.)) 
Hr^ = + «i + nj) + r(4ci + c{c2 + (h.c.)). 

But, from the point of view of the root site, in F2 r is directly 
linked to one site (1), whereas in Fq r is directly linked with 
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FIG. 2: Inequivalent rooted graphs for N — L — 4; the root vertex 
is directly linked with 1 (a), 2 (b) and 3 (c) vertices. 



two sites (1,2); moreover the subgraphs have different con- 
nectivity:^;^,^ = {(1,2)}^^^^ =0. 

The Hamiltonian based on Fi differs from the other two 
since there are three hopping terms. In this case r is directly 
connected to all the other sites as in Tq but here the connec- 
tivity of the subgraph is different, i?p^ = {(1,2)}. 

Among all the possible graphs for the case N = L = 3 we 
do not consider disconnected graphs like the ones in Fig. lb. 
The first case is a particular instance of the system L = 2. In 
the second case, since in our analysis the initial state will al- 
ways be a product state belonging to the site occupation num- 
ber basis, no entanglement can be genertated between the root 
mode and the others. 

Following the same ideas one can derive the inequivalent 
F's also for the case N = L ^ 4. One can start by choosing 
the root site and connect it to just one of the other vertices; 
then one has to find all the subgraphs F' that give rise to in- 



equivalent rooted graphs and repeat the whole procedure for 
the cases in which the root is dkectly connected with two and 
three of the other vertices. The result of this process is given 
in Fig. 2. Once selected the relevant rooted graphs F^ one can 
build the relative Hamiltonians Hr^ (e, r) and the relative evo- 
lution operators (t) — exp{~itHrj (e, r)). The entagling 
power of the latter can then be calculated as described in the 
next section. 



IV. MODE ENTAGLEMENT 

The problem of the entanglement of two-mode (compo- 
nent) Bose Einstein condensates (BEC) has been recently 
studied in different ways |7, 11, 13, 14]. In Ref. iU the 
system is composed by a two component BEC consisting of 
N atoms that can access two different hyperfine states \A) and 
\B) coupled by an electromagnetic field. In this paper it is 
stated that the maximally entangled state for the system is 
given by 

|Gi/Z)A. = i=(|A)®" + |S)«") 

which is the N-partite generalization of the GHZ (Greenberg- 
Horne-Zeilinger) state for 3 qubits 

\GHZ)r, = ^{\000) + \lll)). 

The underlying decomposition of the system is defined by 
considering each of the N bosons as a subsystem living in the 
two dimensional Hilbert space given by span{ | A) ,\B)}. 

In this picture to each boson corresponds a qubit and the 
Hilbert space of the total system then corresponds to the 
Hilbert space of N qubits. On the other hand, as it has 
been pointed out in 1 13], the idea to consider different bosons 
as different two-levels subsystems and to calculate the value 
of the entanglement according to some suitable measure, al- 
thought mathematically possible, it is physically questionable 
in view of the indistinguishability of the particles. In order to 
have some entanglement between two systems, these have to 
be separately accessible from the experimental point of view. 

In general one can think to speak of the entanglement be- 
tween indistinguishable particles, but it is not possible to make 
any use of this resource since there are no physically realiz- 
able measurements on the system that can discriminate be- 
tween the particles. 

In 0] one of the systems under study is a pair of tunnel- 
coupled BEC's considered as a realization of a bimodal BEC. 
In this paper. On the contrary of what is done in Ref ] 14], in 
order to calculate the entanglement present in the ground state 
of the system, the subsystems that are considered physically 
meaningful are the two modes. The number of bosons present 
in the system is always fixed. In this set up it is possible to 
have physical access to the subsystems. For example, one re- 
alization of the model is given by a BEC with a fixed number 
of bosonic particles that are trapped in a two wells potential, 
where each well (site) can be considered as a mode; here the 
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modes are spatially separated. The indistinguishability of the 
particles does not allow to determine experimentally in which 
mode is a specific particle, but it is physically possible to mea- 
sure the occupation number of each mode, i.e. the number of 
particles present at a given instant in the wells. The system 
can be viewed as a bipartite system for which a good measure 
of the entanglement exist: the Von Neumann entropy of the 
reduced density matrix. In this paper we follow this point of 
view and we make use of the notion of mode entanglement to 
characterize our systems. 

We consider systems in which the number of modes (sites) 
L is fixed and is always equal to the number of particles 
N e {3,4}. The Hilbert space ij(^) of the system for a 
fixed = L is a subspace of the Fock space h®^, where 
h — span{|n)}„giN. The subspace can be written as 



i/W ^spanj^jCo^' 



where 



A = {n = (no, ni, .., ni_i) : rij £ IN and 2, "-j — 



L-l 

3=0 



The general state of the system can then be written as ~ 
EsGA"c(ri)|n). 

We are interested in calculating the entanglement of one of 
the modes, say the mode 0, with respect to the rest of the sys- 
tem. In practice we can always reduce to a bipartite picture 
and we make use of the Von Neumann entropy to calculate 
the entanglement of the relevant mode. To do this we need 
to calculate the on-site reduced density matrix. It is interest- 
ing to note that i/^^^ cannot be written in terms of the tensor 
product of the Hilbert space of the mode subsystems: in order 
to calculate the on-site reduced density matrix one has to re- 
mind that the iJ is a subspace of the entire Fock space, that 
by definition is written as a tensor product of the mode sub- 
systems. Thus one can borrow this tensor product structure in 
order to calculate the trace over all the modes but the one. 
Let be V — {0, 1, .., i — 1} the set of sites (modes, or vertices 
in the graphs terminology) then the reduced density matrix of 
the site can be written as: 



P 



(0) 



J2{n\i c{h)c*{k)\h}{k\ 

j=0 \nU) 

N 



where 



A\o = ("•i,--,"-L-i) : e IN}, 
and where — (j, ni, .., e A^. 



In view of the constrain given by the fixed number of parti- 
cles is diagonal and the Von Neumann entropy relative to 
the mode can be simply written as: 



= (0) 



N 

= -E^fiogsP^' 

3=0 



(0) 



That is, e'"^ is the Shannon entropy of the probability distri- 
bution 



(0}^ 



(2) 



In order to characterize the different systems under consid- 
eration, and that are described in detail in the next section , we 
use of the following functional 

EP("^ = [ sup (e(")(t))]/log(Ar+ 1) 

tG[0,T] 

The expression sup^gjo t] (s^'^HO) represents the maximum 
value of the entanglement of site with respect to the rest of 
the system during the evolution of the system; this value has 
been nomalized with respect to log(A^ + 1) that is the max- 
imum value of the mode entanglement that can be reached 
by the system and that corresponds to the situation in which 
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(0) = 



We take the value of EP^"'' as a measure of 
the entagling power which characterizes the pair {Urj , \ijJin))\ 
what we want to compute is, for each Tj, the entangling ca- 
pability of the evolution operator Urj when it is applied for a 
given period of time to the initial state iV'in)- 

As far as the choice of the initial states is concerned, we 
decide to use two particular kind of states for both the case 
N = L = 3 and N = L = 4. The first kind of states are 
the uniformly occupied ones : |111),|1111). The second kind 
of states are the ones in which all the particles are localized 
in the mode 0: |300), |4000). This states allow us to test the 
entangling power of the operators Urj in two opposite and in 
a sense "extremal" situations with respect to the localization 
of the particles in the sites. 

This choice is also motivated by the possible physical im- 
plementations of the model under study. In fact one of these 
possible physical realizations is given by a set of ultracold 
bosons placed in an optical potential. The goal of placing ex- 
actly one atom in each well of an optical lattice created by a set 
of counterpropagating laser beams has been recently achieved 
iHtJ thanks to the realization of a super-fluid/Mott-insulator 
phase transition. This has motivated us to choose, among all 
the possible initial states of the evolution, the ones for which 
the site occupation number is equal to one, {rij = 1, Vj — 
0, .., L — 1). In this kind of application, in order to realize 
the actual systems we are interested in, one should be able to 
isolate the relevant sites with the required graph topology. 

The choice of the states in which all the particles are lo- 
calized in one single mode is motivated by the fact that the 
typical realization of a BBC consists in the localization of all 
the bosons in a trapping (parabolic) potential; with a further 
step it is possible to create a lattice, by means of properly 
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sized counterpropagating laser beams. In this kind of realiza- 
tion one should be able to isolate the relevant sites with the 
required graph topology and to fill with the proper number of 
particles the desired mode (well, site). 

In general, for the realization of our model, one should also 
be able to manipulate the control parameters of the model, 
that is to fix the self-interaction and hopping parameters to the 
desired values, that are the same for all the sites and all the 
links between them. Thus, althought the organization of the 
sites in graphs could be a difficult task, the realization of our 
model requires only a "global" control of the parameters since 
it does not require to perform control operations on single sites 
and/or single links between sites. 



V. RESULTS 

Our goal is to characterize, for each {Urj, iV'm))^ the be- 
haviour of the entangling power EP^^^ (e, r) when the ratio 
between the self-interaction parameter and the hopping pa- 
rameter e/r varies from zero to a value that is much greater 
then one. In order to do so we fix r = 1 (that is we measure 
e in T units) and let s vai-y in [~emax,emax] (Smax > 0) with 
step Ae. For each e £ [—£max, £max] we simulate the full 
quantum evolution for t e [0, T] that the system undegoes ac- 
cording to the unitary operator Upj (t) = exp(— ti/r^ {£, t)), 
{Ti = 1), when the initial state is \ipin)- The evolution time 
T has to be taken much greater than all the eigen-periods of 
the syetems i.e., T ^ muXiE^^ , where Ei denote an energy 
eigenvalue. If this latter condition holds one has no changes 
in the entangling power by allowing a longer evolution time. 
We explicitly verified this asymptotic behaviour in our numer- 
ical simulations. For each of these evolutions we calculate the 
value of EP}!''^ . 

In this way we can consider three different regimes. When 
e = there is no self-interaction between the bosonic particles 
but these are free to hop between the different sites according 
to the graph F. In this regime the intial state |V'm) is not an 
eigenstate of the pure hopping Hamiltonian. When e £max 
we have e/r ^ 1, the hopping between the different sites 
is highly reduced and the systems tend to evolve toward states 
which are very close to the initial ones \ tpin) = |^), that in this 
case are eigenstates of the pure self-interaction Hamiltonians 
corresponding to the eigenvalues e '^fZ^^ nf- The last regime 
is the intermediate one in which there is a competition of the 
two processes that govern the evolution, self-interaction and 
hopping. 

To simulate the evolution of the systems under study we 
use a second order Runge-Kutta routine to integrate the time 
dependent Shrodinger equation for the different Hamiltonians. 



A. N=L=3 

In this subsection we illustrate the results of our simulations 
for the case N = L = 3. In the following we refer to the 
figures 3, 4 and 5. For each graph Tj,j e {0, 1, 2} we have 
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FIG. 3: Plot of £P^g(|e|) for \ipin} = |111> (iandiii) and for 
IV'in) ~ |300) {ii and iv); the curves i and ii refer to |e| = e, while 
the curves Hi and vi refer to |e| — —e. 
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FIG. 4: Plot of EP^^{\e\) for = |111> (i and in) and for 

IV'in) = |300) {ii and iv)\ the curves i and ii refer to [e] = e, while 
the curves Hi and vi refer to |e| = — e. 



computed EPf,^ as a function of e, using Ae = 0.1. For 
each value of e we evolved the system from t = OtoT = 15 
(arbitrary units). For the curves [i) and the [Hi) the intial 
state is \'tpin) = The curves (ii) and the {iv) refer to 

simulations in which the initial state is \ipin) = |300). In the 
{i) and {ii) cases e varies in the interval [0,emaa;], while in 
the {Hi) and {iv) cases e e [—Emax, 0]; for all the different 
cases e„iax — +20 and the value of i?Pp°^ has been plotted 
in function of jej. 
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FIG. 5: Plot of BP^J^del) for = |111> {i and Hi) and for 

I'i/'in) ~ |300) (ii and iv); the curves i and ii refer to |e| = e, while 
the curves in and refer to |e| — —e. 

A first interesting feature about the systems under study is 
that when e = the entagling power of the evolution operator 
does not depend on the graph Tj but only on the initial 

state. In fact, EP'f^ ^ 0.92, Vj, in the case |V^,;„) = |111) 

while EP'^^^ 0.91, Vj, for |?A,„) = |300). The entagling 
capability of the Ur/s , when e = 0, is then only slightly en- 
hanced if the system starts in a situation in which the particles 
are uniformly distributed over the sites. 

We now turn to analize the maximum values reached by 
EPj.'^^ for the different sistems. A first fact that emerges from 
our simulations is that the maximum theoretical value of the 
mode entanglement (log2 4) can be reached for any j if the 
system starts from \-ipin) = |300). In fact, for an appropriate 
choice of the ratio \s\/t, regardless of the sign of the self- 
interaction, 300), e) > 0.99. lf\i^^n) = I 111) this 
is true only in the case j — 0. Another interesting feature, 
common to the cases just mentioned, is that EP^^^ increases 
very rapidly when |e| starts to differ from and passes from 
0.91 — 0.92 to values > 0.97 in correspondance of \e\ — 0.1; 
then the value of EP^^^ remains very high (> 0.95) for 
|e| < e', where e' > differs for the various cases. We 
can the say that, in this interval of values and for the above 
mentioned choice of (Tj, \ijjin)), the self-interaction between 
the particles, either repulsive or attractive, enhances the capa- 
bility of the 's to entangle the mode with the rest of the 
system. 

On the contrary, in the cases (Fi, |111)) and (r2, |111)) 
the maximum theoretical value of the mode entanglement is 
never reached and turning on the self-interaction only slightly 
enhance the entangling capability of the [/r^ s; this happens 
in correspondance of values of |e| which are smaller or of the 



same order of t(= 1). 

Another interesting feature highlighted by our simulations 
is related to the change of sign of the self-interaction parame- 
ter (e — e) and a particular property of the graphs Fj's. If 
we compare the simulations i's vs iii's we can see that for the 
graphs Tq and r2 

EP^fis)=EPl.fi-e) ye e [0,20] 

while for Fi the curves i and Hi differ significantly; the same 
holds for the curves ii and vi. The difference of the behaviour 
of EP^°^ in the two cases depends on a particular property of 
the graphs we have used to characterize our systems. These 
can be distinguished in bi-partite and non bi-partite graphs. A 
graph F is bi-partite if 

F A U S and (a, 5) e 4^ a e A and 6 e S. 

Suppose for example to colour each of the vertices of a graph 
with one of two different colours A and B. A graph is bi- 
partite if the vertices coloured with A are linked only with 
vertices coloured with B, and viceversa. In our case the bi- 
partite graphs are Fq and F2 (note that the bi-partiteness of 
a graph depends only on the pair (V, E) and is independent 
from the root). 

If Tj is bi-partite, the Hamiltonian that results from the 
change of sign of e can be written as 

i/r,(-£,T) = -P-ifr,(£,T) -Pt, (3) 

where P = cxp(—i7r ^^^^ rij ) is the unitary operator that 
implements the canonical transformation 

^ (-l)'^-"^^'^Cj (4) 

(XA is the characteristic function on the subgraph A). The 
evolution operator for a given bi-partite graph can then be ex- 
pressed as 

C/r/(t) = cxp(-i<i/r,(-e,T)) =Pexp(iti/r,(e,T))pt 
= PU^^i-t)P^ 

If we now apply this operator to a state belonging to the 
occupation number basis \n) we obtain: 

iV'rW) = t/f/WI"> =-Pexp(iti/r,(£,T))exp(~i0«)|n) 
= ^ exp(i0^)as,7i(-i)exp(-i0s)|m) 

mGA" 

where exp(— i^fj) and exp{i(f)^) are, respectively, the 
phase factors given by the action of the operator P^ on the 
initial state \n) and of the operator P on the states |m) in the 
sum. The terms t) can be obtained by acting with the 

time reversed operator Uf-. {—t) on P^|n) when it is written in 
terms of the energy basis states, in which the Hamiltonian is 
diagonal, and then going back to the site occupation number 
basis. Since Vj Hpj is a real simmetric matrix there exist a 
real orthogonal matrix O (O ■ = 1) that diagonalize it and 
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with which it is possible to implement the required changes of 
basis. We can then write 



k 



Pnmkexp{+itEk), yn,m 



where Ek are the eigenvalues of the Hamiltonian and (3nmk 
can be written in terms of the real entries of O. This fact is 
crucial when we compare the evolutions of \n) obtained by 
applying U^^{t) and Ul..{t) (note the sign of t). The latter 
gives: 

\rn{t)) = eM-^tHr^ie,T))\n) 



where 



Q-nm{t) =^Pfimkexp{-itEk), \/n,rn. 



We have now the elements to explain why, when Tj is bipar- 
tite, the curves {i) and (ii) overlap perfectly with the curves 
in) and vi) respectively. Since the terms Pftmk are real we 
have that for any t and Vn, m the amplitudes of the states |m) 
in the expression of \tjjZ'^{t)) and |V'|(*)) satisfy: 

I exp{i(j)^)afiA{-t) exp{-i(j)fi)\'^ = |a^^A(^)|^ (5) 

in fact a,i^{~t) = ce^,fi{t). Consequently, since the value 
of the mode entanglement and of EP^^ depends only on 
the modulus square of the probability amplitudes of the ba- 
sis states |m) in the expression of \^'^^[t)) (see section IV), 
for any given ji/im) — \n) and for all the bi-partite graphs we 

have EPf^{£) = Ep!f^ {—e). In general, this is no longer 
true for the non-bipartite graphs. 

Another relevant and general feature of the entagling power 
VP's and for both the selected initial state is that EP^^ rapidly 
decreses when |e| is significantly greater than 1(= r) with 
a rate that is higer in the case |V'in) = |300). This can be 
explained with the following argument. When |e| becomes 
significantly greater then 1, the Hamiltonian of the system, 
V Tj, can be considered as composed of an unperturbed diag- 
onal term H^, the self-interaction one, and a perturbation term 
given by the hopping matrix Hr- Both the initial states of our 
computation are eigenvectors of the unperturbed term and cor- 
respond to the eigenvalues 3e (|111)) and Qe (|300)). As |£| 
increses, the perturbation becomes smaller and the only pro- 
cesses relevant for the evolution are the first order ones. This 
means that the probability distribution [p^^ , pf^ , p^^ , p'f^ ), 
given by the diagonal elements of is higly unbalanced and 
this results in a value of e(°^(i) that for any t is very smaU, 

and so it is, consequently, the value of -EPp°^ . 

To give a more detailed explanation of this phenomenon 
let's take for example To and suppose |V'm) = I HI)- Iri 
general, at first order, the state |V'm) can be coupled by the 
perturbation to the states |n') S v where 



{|n') I n'eA^and(n'|iJ,|^i„)7^0}. 



In our case x\iii) = {|021), |012), |201), |210)}: the modu- 
lus square of the probabihty amplitudes of these states con- 
tribute to the terms p^^^ , p^^ of the distribution. The state 
in which no = 3 is coupled to |111) only by second order 

processes and for this reason the term ^3°^ <C Po°^ , Pi°^ , ^2°'' ; 
furthermore, when the perturbation is very small, the terms 
Po°^ , p'^^ <C pi'^ and the distribution is very peaked around 

(0) 

Pi ■ 

The same arguments hold for Fq when = |300), but 
there is a significant difference between the two cases. In fact 
X|300> ~ {|201), |210)} so that when the pertubation is small 

Po°^ , p^i^ -C P2°^ <C p^^ : the probability distribution is more 
unbalanced than in the previous case and this gives a first ex- 
planation of the fact that in the perturbative regime 



i;p(f(|iii))>EP^^^(|3oo)), vj. 

There is another important difference between the two cases 
that justifies this relation: the strenght of the coupling. This, 
at first order, depends on the ratio r/Ai?, where Ai? is the 
difference in energy between IV'm) and the states belonging 
to Xj^. )• The ratio is higher for I -(/jin) = (l/2|e|)then 
for \ipin) = |300) (1/4|£|) and this results in higer values 
of probabiUty relative to the states in xj^ ) in the first case. 

Consequently the associated eigenvalues of p^"^ can be higher 
for IV'in) = I 111) thus giving rise to a probability distribution 
which is flatter and to higher values of EPj 



(0), 



(0) 



We now use the perturbative approach to describe the fea- 
tures relative to the non-bipartite graph Fi. Here the inter- 
esting thing to notice is again the role of the sign of the self- 
interaction, which depends on the initial state. For |i/'m) = 
[300) a repulsive (positive) self-interaction enhances the en- 
tanghng capability of Uvj {t). In fact, for a large part of the 
interval [0, |£max|], 



EP 



(0) 



(|300),+|£|)>i?p(°)(|300),-|£|) 



(6) 



On the other hand, when |i/;i„) — |111) is the attractive (neg- 
ative) self-interaction that favors the entangling power of the 
evolution operator. In fact for any |£| € [0.5, Smax] 



EP^y\lll),-\e\)> EP^> {\lll),+\e 



(0), 



(7) 



This behaviour can be understood in detail at least in a per- 
turbative picture. Let's start from the case |^i„) = |111). 
When |£| is significanly greater than 1 (= r) we can use a two- 
level picture to describe our system. The first level is given by 
\i)in), while the second one is given a uniform superposition 
of the g = Ixj^j^) I degenerate states that are coupled to |'^i„) 
by first order processes, that is 

|J) = (1/V9) E \n')^Hr\^in). (8) 



The two-levels Hamiltonian can be written in general as 



H : 



El W 
W E2 



(9) 
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where 

El = {ll>in\He\tpin) , 

W = {I\Hr\i'^n) = {An\Hr\I)*, 

E2 = {I\H,\I) + {I\Hr\I) . 

It is a standard result of two coupled levels systems theory that 
the probability of finding the system in the state |/) at time t 
can be written as (Rabi's formula) 

P\^,^)^\i){t) = |(7|expHt(ffe + ifr)|V'm)|^ 
A\W\^ 



4|W^|2 + {Ex - E2Y 



where 9 = a/4|M^|2 + {Ei - £^2)^/2. Here the important 
fact is the role of the term [Ex — E2Y in the denomina- 
tor from which depends the maximum value attainable by 
P|^.^)^|/) {t)\ this is directly Unked with the maximum val- 
ues of p'^^ , and consequently with the maximum value 

attainable by EP!^^ in the perturbative regime. In fact it is 
possible to reproduce, with a good approximation, the results 
of our simulations by using this probability. This can be done 
by calculating the mode entanglement of the following state: 



\/r^iV'm) + Vp|/) 



where 



P{e,T) 



4|Ty|^ 



A\W\-^ + {Ex - E2f 



(10) 



(11) 



is the maximum attainable value of P\^,.^^^\j^ {t) according to 
the Rabi's formula. In the case of the system Fq, the state |/) 
is given by a superposition of states for which no = 0, 2. The 
only contribution to the terms p"^^ , p^p of the reduced density 
matrix of the mode are given by these states. This means, as 
described in the previous paragraphs, that the higher are the 
values of these terms the "flatter" is the probability distribu- 
tion (pq"^ , , ^2°^ , p^^'' , ) and consequently the higher is the 
value ofi;P]l°^Forthis reason the state \ (b) can be considered 
as an example of state in which the entanglement is maximum, 
and we can compute the entangling power for any rooted 
graph, and for all the values of e in the perturbative regime 
as£P^°^ =-Tr(pW log2p(°)) where = Try^J|.;i)(^|). 

By means of the two levels picture we now have the oppor- 
tunity to see from another perspective the difference between 
bipartite and non-bipartite graphs. This result is essential to 
explain the difference in the behaviour of EP^^ described by 
equations (6) and (7). Since both \iljin) and |/) are eigenvec- 
tors of He, if {I\Hr\I) = 0, the difference {Ei - E2) does 
not depend on r and is linear in e; thus the sign of the lat- 
ter does not have any effect on P. This is always true for 
all bipartite graphs. To demonstrate that for all bipartite graphs 
{I\Ht\I) = 0, we first observe that, since |/) oc HT\tpin), we 
can always write 



We can then make use of the transformation P = 
exp{—TTj2j£A '^j) previously described to demonstrate that 
for any bipartite graph (lll|i?^|lll)=0. In order to do this 
we calculate the expectation value ^lll|Pti/^P|lll^ in two 
equivalent ways; in the first we let act P(P^) on |111)((111|) 
to obtain: 

((lll|pt)i/3(p|lll)) = (lll|H3|lll). 

In the second one we use the fact that for any bipartite graph 
P^HrP = -Hr, this allows to write 

(iii|pi'i?^p|iii) = _ {in\H^\ni) . 

This imphes that for any bipartite graph 

(/|l/,|J) = (lll|if3|iii) = o. (13) 

The last result is, in general, no longer true when the graph 
is non-bipartite. In fact, for Fi, {I\Ht\I) is constant and pro- 
portional to r. This is essential to explain the difference in the 
behaviour of PP^"^ in correspondance of an attractive or a re- 
pulsive self-interaction. In fact, the presence in E2 of the term 
proportional to r renders the difference {Ex — E2) sentive to 
the sign of e. 



(^i-^2) = -(2e + 3r). 



This implies that 



P{-\s\)>P{+\s\), 



and consequently, as one can see by calculating EP^^ from 
the state \6), 



Pp(°)(|lll),-|e|)>PP^;'^(|lll),+|e|). 

Since the resuts (12) and its consequence (13) hold also 
for \ipin) = 1 300) we can use the same arguments to give a 
justification of (1). Here simple calculations give: 

Ex - E2='ie- T. 

The role of e is then opposite with respect to the previous case: 

P{+\e\)>P{-\e\), 
and consequently, in the perturbative limit 

EP^y\300),+\s\) > i;P^°^(|300),-|£|). 
B. N=L=4 

In this subsection we illustrate the results of our simulations 
for the case N = L = 4. In the following we refer to the 
figures 6-16. For each graph Fj, j € {3,4, .., 13} we have 

computed EPjP^ as a function of e, using Ae = 0.1. For 



(0), 



{I\Hr\I) OC (V'inl-ffxIV'm) = • (12) each value of e we evolved the system from < = to P 
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FIG. 6: Plot of EPr^ (|e|) for \ipin) = |1111) (i and m) and for 
I'l/'in) ~ |4000) (ii and ill); the curves i and ij refer to |e| — e, 
while the curves in and vi refer to |e| — —e. 



FIG. 8: Plot of EPr^{\e\) for |Vin> = |1111) [i and m) and for 
IV'jn) = |4000) (iiandin); the curves i and ij refer to |e| = e, 
while the curves Hi and vi refer to [e| = —e. 




FIG. 7: Plot of _BP^^ (|e|) for |Vin) = |1111) (i and Hi) and for 
IV"!") = |4000) (iiand in); the curves i and ij refer to |e| = e, 
while the curves Hi and vi refer to |e| = — e. 



FIG. 9: Plot of EP^^{\e\) for I^Z-in} = (i and m) and for 

IV'in) = |4000) (ii and ill); the curves i and ii refer to [e| = e, 
while the curves Hi and vi refer to |e| = —e. 



15 (arbitrary units). The curves (i) and {Hi) in each graph 
are the resuhs of the simulations for \%l]in) = |1111). The 
curves {ii) and {iv) refer to the simulations in which the initial 
state is IV'm) = |4000). In the [i) and [ii) curves e varies 
in the interval [Ojemaa;], while in the [Hi) and [iv) curves 
e E [—emax, 0]; for all the different cases Emax = +20 and 
the value of EP^^^ has been plotted in function of |e|. 

As we did in the previous section we start by observing 
that when e = the entagling power of the evolution oper- 



ator Ur^ does not depend on the graph Tj but only on the 
chosen initial state. In fact i?Pp"' = 0.83 , for any j, in the 
case iV'in) = |1111), while EP^°^ = 0.87 , for any j, when 
iV'in) — 1 4000). We see here a first difference with case 
N = L — 3: the entagling capability of the Ur/s , when 
e = 0, is enhanced if the system starts with all particles local- 
ized in the mode 0. 

We now turn to analize the maximum values reached by 
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FIG. 10: Plotof £P^j(|e|) for \ipin) = |1111) (i and ui) and for 
IV'in) = |4000) (ii and ill); the curves i and ij refer to |e| = e, 
while the curves Hi and vi refer to |e| = —e. 




5 10 15 20 



FIG. 11: Plotof BP^jdel) for = |1111) (i and iii) and for 

I'i/'in) = |4000) (ijandin); the curves i and ij refer to \e\ = e, 
while the curves Hi and vi refer to |e| = ~e. 



EPr for the different sistems. The maximum theoretical 
value of the mode entanglement (log2 5) is always reached 
V Tj when the initial state is |4000). In fact, for an appropriate 
choice of the ratio |e|/T, regardless of the sign of the self- 
interaction, EPjP^lAOOO)) > 0.99 . In contrast with the case 
N — L — 3, none of the systems reaches this value when the 
initial state is the uniformly occupied one; the highest values 
of EP^°\\1111), ) are > 0.97 and they are reached by the 
systems Fs, Fn ( only for e > 0), and Fiq. 




»-» 1 TT-rm- i 

5 10 15 20 



FIG. 12: Plotof £P^^(|e|) for \1p^n) = |1111) {i and Hi) and for 
IV'jn) ~ |4000) {ii andiv); the curves i and ij refer to |e| — e, 
while the curves Hi and vi refer to [e| — —e. 




FIG. 13: Plot of EP^^ (|e|) for \ip,„) = [HU) (i and Hi) and for 
lil^in) ~ |4000) (ii and ill); the curves i and ii refer to je] = e, 
while the curves Hi and vi refer to |e| — — e. 



When l^in) = 1 4000), EP^^ ,' increases when |e| starts to 
differ from and passes from 0.87 to values > 0.92 in cor- 
respondance of \e\ = 0.1; then the value of EPjP^ remains 
very high for \e\ < e', where e' > differs for the various 
cases. We can then say that, in this interval of values the self- 
interaction between the particles, either repulsive or attractive, 
enhances the capability of the Urj 's to entangle the mode 
with the rest of the system. 

In the case IV^m) = as mentioned above, the max- 
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FIG. 14: Plotof £P^j(|e|) for \ipin) = (i and ui) and for 

I'l/'in) = |4000) (ii and ill); the curves i and ij refer to |e| = e, 
while the curves Hi and vi refer to |e| = —e. 



FIG. 16: Plot of £P^^(|e|) for |?/>»„) = |1111) (i and ui) and for 
IV'jn) = |4000) (iiandin); the curves i and ij refer to |e| = e, 
while the curves Hi and refer to [e| = —e. 




FIG. 15: Plotof £;P^\(le|) for = |1111) (i and wi) and for 

IV'ii) = |4000) (ii and ill); the curves i and ij refer to |£| = e, 
while the curves Hi and vi refer to |e| = — e. 



imum theoretical value of the mode entanglement is never 
reached and the effect of turning on the self-interaction de- 
pends on the system and on the sign of the self-interaction. 
However, the enhancement of the entangling capability of the 
's, when present, is found in correspondance of values of 
\e\ which are smaller or of the same order of r(= 1). 

By observing the figures 6-16 we can see that also for the 
case iV = L = 4 we can distinguish between two different 
behaviours of Ep!^^ with respect to the sign change of e. For 



j e {3, 4, 6, 8, 10} we have that for both the intial states 

EP^°\e) = EP^°\~e) Ve e [0, 20], (14) 

that is the curves i and ii perfectly overlap with the curves 
Hi and iv respectively. On the other hand for j E 
{5, 7, 9, 11, 12, 13} this is no longer true. 

This fact can be explained using the same arguments we 
used in the previous section. In fact the Tj,j — 3, ..13, can be 
distinguished in bi-partite and non-bipartite graphs. For the 
bipartite ones we can apply to the related Hamiltonians the 
unitary transformation (4) in order to get equation (3). The 
latter and the subsequent reasonings hold for any N — L and 
thus equation (5) implies that Vj G {3, 4, 6, 8, 10}, i.e. for all 
bipartite graphs equation (14) must hold. 

We now turn to examine the behaviour of EP^^ in the per- 
tubative regime. We first observe that for both the initial states 
EP^^ rapidly decreses when jej is significantly greater than 
1(= t). This happens for the same reasons we described for 
the case N = L = 3. The Hamiltonian, VFj, is composed of 
an unperturbed diagonal term H^, the self-interaction one, and 
a perturbation term given by the hopping matrix Ht- Both the 
initial states of our computation are eigenvectors of the unper- 
turbed term and coiTespond to the eigenvalues 4£ (| 1111)) and 
16e (|4000)). As |e| increses, the perturbation becomes small 
and the only processes that become relevant for the evolution 
are the first order ones. At first order, the state IV'm) can be 
coupled by the perturbation to the states \n') S x\,p- )■ This 

implies that for j-fAin) — we have that p'"^\ p^^'' <C 

Po\p2^ < Pi^^ while for iV'jn) = |4000) we have that 
Po^ , pf'' , p'"2^ <^ P3°^ <^ p'^^^ . The probability distribution 
(Pq*^' , p^°^ , pj""* 1 P3°^ 1 pf^ ) is higly unbalanced in both cases. 
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and this results in the decreasing behaviour of EP^ ,', but it is 
more unbalanced in the case \ipin) = |4000) , Vj, and this is a 
first explanation of the fact that in the perturbative regime 

EP^°^{\1111)) > EP{.f{\mO)), Vj. (15) 

A further explanation for the last relation is given, just as for 
the case TV = L = 3, by the fact that the strenght of the 
first order coupling between |i/'m) and the states belonging 
to x'\^- ) is higher in the case \4>in) — ('''/Ai? — 

l/(2|e|), then in the case l^/'in) = |4000), (r/AS = 1/(6|£|). 
As explained in the previous section this affects the highest 
reachable value for the probability related to the states belong- 
ing to in X| ^ ^ and consequently affects the values reached by 

EP^°^ in the way described by (15). 

We finally come to the description of the role of the sign 
of the self-interaction when |e| is significantly greater than 
1(= r) in the case of the non-bipartite graphs Tj, j e 
{5, 7, 9, 11, 12, 13}. The phenomenology is analogous to the 
one found for the case N = L = 3. For \ipin) = |4000) 
a repulsive (positive) self-interaction enhances the entangling 
capability of Urj{t), that is for a large part of the interval 

[0, \emax\] 

£;F]1°^(|4000),+|£|) > £;P]^°^ ( 14000), (16) 

In contrast, for jf/'m) = IHH) it is the attractive (negative) 
self-interaction that for a large part of the interval [0, |ema2;|] 
enhances the entangling power of Urj (t): 

£;p(°)(|llll),-|e|) >i?p(f (11111), +|e|). (17) 

This behaviour can be understood by means of the same argu- 
ments exposed for the case N — L — 3 that we briefly recall 
here. When |e| is significantly greater than 1(= t), we can 
again use a two-levels picture to describe our systems. The 
first level is |V'm), while the second |/) is given by a uniform 
superposition of the q — |x'^ ^| degenerate states that are 
coupled to \^in) by first order processes, see eq. (8). The 
two-levels Hamiltonian is again given by (see formula) and 
the probability P\^>i„)^\i){t) is expressed by the Rabi's for- 
mula. As we have seen for the case = L = 3 with the 
introduction of the state |0), the maximum attainable value of 
this probability P = P{e, r) is directly related to the value 
of EPj?\ that is the higher the former the higher the latter. 
In this context it is again the term {Ei ~ E2) that is crucial 
to give an explanation of the equations (16) and (17). In fact, 
equation (13) holds also for the case N = L = A and for the 
intial states |1111) and |4000). Thus, for the bipartite graphs, 
wehavethat (/|i/T-|/) = 0, the term (i?i — P2) is independent 
of T and depends linearly on e. This implies that the proba- 
bility P|^,.^^^|/^ [t) does not depend on the sign of e for the 
bipartite graphs Tj, j e {3, 4, 6, 8, 10}. 

In the case of non-bipartite graphs this is no longer true and 
the dependance of (Pi — P2) on r give rise to a sensitivity to 
the sign of e of P = P(e, t). Hence the different behaviour 
of EP^^ described by equations (16) and (17). 
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FIG. 17: Plotof £Pr,(|e|) for j = 5,7,11; |V'»n) = |1111) and 
|e| = e. 
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FIG. 18: Plotof £;P^^.(|e|) for j = 10,11,12,13; = |1111) 
and |e| = e. 



C. Graph Properties and Entangling Power in the Pertubative 
Regime 

In the last part of the paper we want to deepen our analy- 
sis in the perturbative regime in order to describe the existing 
relations among the behaviour of EPjP^ for the different Tj. 
As we have seen in the previous paragraphs, in the perturba- 
tive regime, this behaviour can be described by using for our 
systems a two levels approach. In this picture, the maximum 
value attainable for PPp"' is directly linked with P, i.e. the 
maximum of the probability P|,0.^^^|/^ (t), see equations (8)- 
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i^mCn such that 




FIG. 1 9 : Plot ofEP^^(\e\) for j 
and lei — — e. 



10,11,12,13; |V'»n) 



mil) 



(11). It is possible to reproduce, with a good approximation, 
the results of our simulations by using this probability and by 
calculating the mode entanglement of the state (10), that is: 



£;p(°)=-Tr(p(°)log,p(°)) 



(18) 



where p*^*^) = Trvy„{\4>) (</>!)■ In order to evaluate (18) we can 
use analytical expressions for P and for p^*'^ that allow us to 
Unk EpI°^ to some of the topological properties of the graphs. 
These expressions can be derived for the general case N = L. 
We first start our analysis by considering \ipin) — 



in this case the expression for P is: 



P = 



16t2A: + (2e + t18;3/2A:)2 ' 



(19) 



where IGr'^k = 4|T-y2g| = 4|M^|2 is given by the cou- 
pling induced by the hopping term Hr and is proportional 
to fc = \E\ = q/2 = \x\^,. ^1/2, i.e. the number of links 
of the graph; 2£ = {I\Hs\I) - {ipin\He\ipin)- The term 
T^lSls/q — T^18l3/2k = {I\Hr\I) depends on fc and also on 
another property of the graph which is the number ^3 of closed 
loops of length 3 present in the graph (one of such loops can 
be seen for example in Fn where the loop is given by the links 
between the sites 0, 1 and 2). 

We now give a brief explanation of the appearance of the 
number ^3 in the expression in the expression of P. As we 
have already pointed out in section V.A, eq. (12), {I\Hr\I) 
is proportional to (?/;j„|i/^|?/'i„). The fact that this term dif- 
fers from zero only if the graph is non bi-partite can be now 
explained in terms of the presence of loops of length 3 in the 
graph. In fact a graph is non bi-partite if and only if it con- 
tains at least one of these loops. The presence in the graph of 
one of such loops results in the presence in iJ^ of terms like 



{'^in\'-^i,j.h,k,m,n\'^in} 0- 



(20) 



For example, in the case Fn, CqC2c\cic\cq is one of such 
terms. In order to calculate {I\Hr\I) we first have to com- 
pute the number of terms for which equation (20) holds. This 
can be done, for any loop of length three, as follows. We first 
observe that the number of possible choices for the index n is 
3, that is the number of vertices involved in the loop. Since 
each of these is linked to the other 2 vertices in the loop we 
have 2 possible choices for to. Once n and to are fixed there 
are only two kind of terms that satisfy equation (20), that is 
a) Tn,h.h,va.m,n and b) Th,m,n,h,ni,n, whcrc h is the third re- 
maining vertex. Consequently, for any loop of length three we 
have 12 terms that satisfy equation (20). By inspection we see 
that for the six terms of type a) {tpin\Tn,h.h,ni.m,n\'^zn) = 4, 
while for the other six {i'in\Th,m,n,h,ni,n\i'in) = 2. This 
implies that, taking in account of the fixed contribute given 
by each of the loops of length three present in the graph, 
{I\Hr\I) = m/q. 

We have then succeded in linking P to the graph properties k 
and ^3.The next step is to explicitly write i?Pp"^ as the entan- 
glement of the state |(/)); to do this it is necessary to express 
Pq"-* , pf^ and in terms of the probability P. This can be 
done in a way that takes in account of the properties of the 
graph Tj and of the state |/): 



Pi 



(0) 



(l-P) + -2(fc-fco) = l-^- 



„(o) (0) 



p. 



Pko 



(21) 



(22) 



where q — ^| — 2k; fc is the number of links in the 

graph and fco is the number of links starting from the root ver- 
tex. The term ■^2(fc — fco) in (21) takes in account of the 
2(fc — fcg) states in |/) that are coupled by first order pro- 
cesses to IV'm) and for which uq = 1. We can now use the 
formulas just derived to describe some of the relations exist- 
ing among EPjP^ for the different rooted graphs in the case 

We first focus on the relation among the values of EPj?\ 
in the perturbative regime, for rooted graphs that differ only 
for the choice of the root vertex. As an example we take the 
rooted graphs F5, F7 and Fn; in figure (17) we have plotted 
the results of our simulations for the three graphs in the case 
iV'in) = 1 1111) and e > 0. From this figure we can see that 



EPj°'^ > EPj. 
ill i 



(0) 



> EP 



(0) 



That is, the higher the value of fcg, i.e. of the coordination 
number of the root vertex of F^ , the higher the value of EPjP '' . 
This is an example of a feature that holds in general for any 
value of = L. In fact, suppose to have 
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• a set of rooted graph A = {Tj } that differ one from the 
other only for the choice of the root vertex; 

• IV-m) = ®jCo^|i),-; 

• a fixed value of £ (> or < 0) in the perturbative 

regime, 

then the Fj's for which _EPp°^ is higher are those for which 
fcg is maximum. Furthermore it is possible to say that 



are and pj^Li- These can be simply written as: 



if > fcj then EP^°^ > EP^ 



5(0) 



(23) 



This can be seen from equations (19), (21), (22). In fact, for a 
given value of e in the perturbative regime, P is the same for 
all the rooted graphs in A, since for all them k and ^3 are the 
same. On the contrary, the flatness of the probability distribu- 
tion (2) and consequently the value of EPjP^ depend, in the 

way described above, on the value of k^. 

We now consider a different situation in which it is mean- 
ingful to compare the values of EPj.^^ for different rooted 
graphs. We take as example the set of graphs B = 
{rio,rii,ri2,ri3}. it is interesting to notice that the non- 
bipartite graphs rii,ri2, and Fia can be seen as generated 
from the bipartite graph Fio by adding respectively one, two 
and three links, that is by increasing the connectivity of the 
subgraph £". It is then interesting to see what happens to the 
maximum attainable value of the mode entanglement when 
the connectivity of the subgraph changes. We refer to the fig- 
ures (18) and (19), which describe the cases e > and e < 
respectively (|f/'m) = |1111))- In figure (18) we can see that, 
for a given £ > 0, the increase of the connectivity of the sub- 
graph E' has the effect of reducing the value of i?Pp°^ . On 
the contrary, when £ < the effect is opposite, see figure (19). 
What one could conclude by looking at equations (21)-(22) is 
that EP^^^ should decrease when the connectivity of the sub- 
graph increses, regardless of the sign of the self-interaction, 
in view of the dependance of pj"^ and p2°'' on ko/k. But this 
contradicts the behaviour described in figure (19). In fact, it 
is the role of the sign of the self-interaction in P that allows 
iJPp°^ to increase with the increase of the connectivity of the 
sub graph E'. We can see here the role of the interplay be- 
tween the self-interaction and the tunnelling processes. 

We can now end our analysis by considering the case N = 
L and l-tpm) = \N) 'S>f^Q \0)j. As a first step we can write the 
analitical expression of P in terms of the graph properties as 
follows: 



P = 



Ar^koN 



Ar^koN + [2{N ■ 



(24) 



where ko is the number of links starting from the root vertex, 
while Z3 is the number of loops of length three that include 
the root vertex. Since for all the states in 1 7) no = N — 1, 
for the determination of EP^^^ in the perturbative regime the 
only terms of the probability distribution (2) that are relevant 



(0) 

Pn 

(0) 

Pn-1 



1-P 
P. 



With the help of (24) we can now see that the result (23) hold 

also for the state \tpin) = \N) 'S'f^Q \0)j- In fact, here it is de- 
pendence of P on ko that determines the behaviour described 
by (23). Furthermor, also in this case it is possible to say that 
the Tj's for which EP^^^ is higher are those for which fc^ is 
maximum. 

Finally we consider again the case B = 
{Fio, Fii, ri2, Fia}, i.e. the rooted graphs that can be 
can be seen as generated from the bipartite Fio by increasing 
the connectivity of the subgraph E'. In this case the increase 
of the connectivity of the subgraph plays a role only if it 
implies the presence in the graph of loops of length three 
that include the root vertex. For the set B we have that 
?3 = 0, 1,2,3 for Fio,rii,ri2, and ri3 respectively. The 
interplay between the self-interaction and the tunnelUng 
processes plays a role which is opposite with respect to the 
case \tpin) = In fact, in the case \ipin) = |4000) we 

have that for a given £ > 0, the increase of the connectivity 
of the subgraph E' has the effect of increasing the value of 
i?Pp°' . On the contrary, when £ < the effect is opposite. 



VI. CONCLUSIONS 

We have analyzed by means of numerical simulations the 
mode-entanglement properties of of bosonic particles hopping 
over graph structures. Starting from an initial unentangled 
state bi-partite mode-entanglement of a graph vertex with re- 
spect the rest of the graph is generated, by allowing for a finite 
time a tunneling along the graph edges. We characterized the 
entangling power of a given hopping hamiltonian by the max- 
imum achieved during the dynamical evolution by this entan- 
glement. We studied the dependence of this entangling power 
by the local self-interaction parameters for all the graphs up 
to four vertices and for two different natural choices of the 
initial state. Our results show the important role of that graph 
topology has in optimization of entanglement production 

The main steps and results,for both cases -/V = L = 3, 4 , 
can be sunnmarized as follows. 

• We have classified the way in which a mode (bosonic 
site) can be linked to the "rest of the world" (other 
bosonic sites) in terms of rooted graphs which are in- 
equivalent from the poin of view of the mode entangle- 
ment of the reference mode (root). 

• When the self-interaction term is absent (= 0) the value 
of EP^^^ does not depend on the graph choosen but 
only on the initial state. 

• We have established for which of the configurations 
(Tj, \ipin)) the free evolution of the system can lead to 



the optimal value of EPj. 



(0) 
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• We have given a general (analytical) justification of the 
different behaviour of EP^^ with respect to the change 
of sign in the selfinteraction parameter. This has been 
done by first observing that the graphs Tj can be classi- 
fied in bi-partite and non bi-partite. We have thus given 
an analitical justification of the fact that 

EP^°^{e) = EPj^fi-e) Ve e [0,20] 

for all the bi-partite graphs. This relation, as it is high- 
Uted by our simulations, is not true for the non bi-partite 
ones. 

• We have given a general explanation for the reduction 
of EPjP'^ when the self-interaction term becomes suf- 
ficiently greater than the hopping term. This has been 
done by adopting a two-levels system perturbative ap- 
proach. 

• For the non bi-partite graphs, in the perturbative two- 
levels picture, we have given a detailed description of 
the following interesting features highlighted by our 
simulations. In the perturbative regime, the entangling 
power of the pair {Urj , IV'in)) is: enhanced (decreased) 
by a positive (negative) self-interaction when, in the ini- 
tial state, the particles are all localized in the reference 



mode (no = N); the opposite holds true when the initial 
state is uniformly occupied {rii — 1, \/ i ~ 1, ..L). 

• In the context of the perturbative two-levels picture 
we have derived analytical formulas that allow to 
link EPjP^ to the topological properties of the rooted 
graphs; this has been done in the general case N ^ L 
and for both the initial states |^i„) = |A^)®^Jo^ |0)j and 

\'>pin) — 'Xi^^q'^II)^. Starting from the simulations and 
with the help of this formulas we have compared and 
described the behaviour of EPj?^ for different graphs. 

We believe that a natural -though not unique - experimental 
framework in which the theoretical analysis carried over in 
this paper can find application is provided by systems of cold 
bosonic atoms in optical lattices. 



VII. ACKNOWLEDGMENTS. 

We are grateful to Radu lonicioiu for useful comments. 
Special thanks are due to Vittorio Penna for introducing us 
to the problem of bosonic lattices and useful discussions. We 
gratefully acknowledge financial support by Cambridge-MIT 
Institute Limited and by the European Union project TOPQIP 



[1] See, e.g., D.P. DiVincenzo and C. Bennet, Nature 404, 247 

(2000) and references therein. 
[2] P. Zanardi, Ch. Zalka, L. Faoro Phys. Rev. A 62 030301(R) 

(2000) 

[3] G. Vidal, I. Cirac, N. Linden, S. Popescu , Phys. Rev. Lett. 87, 
137901 (2001) 

[4] M. S. Leifer, L. Henderson, N. Linden, Phys. Rev. A 67, 012306 
(2003) 

[5] M. Nielsen et al quant-ph/0208077l 

[6] Chris D. Godsil, Gordon P. Royle, Algebraic Graph Theory, 

Graduate Texts in Mathematics, Springer Verlag (2001) 
[7] G. J. Milburn et ai, Phys. Rev. A 55, 4318 (1997); 
[8] D. Jaksch et ai, Phys. Rev. Lett. 81, 3108 (1998). 
[9] Z.B. Chen and YD. Zhang, Phys. Rev. A 65, 022318 (2002) 
[10] R. lonicioiu and P Zanardi, Phys. Rev. A 66, 050301 (2002) 
[11] L.-M. Duan, J. I. Cirac, and P ZoUer, Phys. Rev. A 65, 033619 
(2002) 

[12] C. Simon, Phys. Rev. A 66, 052323 (2002) 

[13] A. P Hines et al. , Phys. Rev. A 67, 013609 (2003). 

[14] A. Micheli, D. Jaksch, J. I. Cirac, and P Zoller, Phys. Rev. A 

67, 013607 (2003) 
[15] L. You, Phys. Rev. Lett. 90, 030402 (2003) 
[16] L.-M. Duan, E. Demler, M. D. Lukin, cond- mat/02105641 



[17] M. Greiner et al. Nature 415, 39 (2002). 

[18] P Zanardi, Phys. Rev. Lett. 87, 077901 (2001) 

[19] P Zanardi, Phys. Rev. A 65, 042101 (2002); P Zanardi, X-G. 

Wang, J. Phys. A:Math. Gen., 35, 7947 (2002) 
[20] J. R. Gittings and A. J. Fisher, Phys. Rev. A 66, 032305 (2002) 
[21] S. van Enk.Phys. Rev. A 67, 022303 (2003) 
[22] Phys. Rev. A 67, 024301 (2003) 
[23] V. Vedral, quant-ph/0302040 

[24] J. Schliemann, D. Loss, and A. H. MacDonald, Phys. Rev. B 63, 
085311 (2001); J. Schhemann, J. I. Cirac, M. Kus, M. Lewen- 
stein, and D. Loss, Phys. Rev. A 64, 022303 (2001). 

[25] R. Paskauskas, L. You, Phys. Rev. A 64, 043210 (2001) 

[26] Y. S. Li, B. Zeng, X. S. Liu and G. L. Long, Phys. Rev. A 64, 
054302 (2001) 

[27] R. Burioni et al, Europhys. Lett. 52, 251 (2000) 

[28] P Zanardi, quant-ph/0212120 

[29] M. P A. Fisher et al., Phys. Rev. B 40, 546 (1989). 



